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Abstract 

(N . 

In this paper, under the hypothesis that p is upper bounded, we construct a Lyapunov functional 
for the multidimensional isentropic compressible Navier-Stokes equations and show that the weak 
solutions decay exponentially to the equilibrium state in L 2 norm. This can be regarded as a gen- 
eralization of Matsumura and Nishida's results in |23| . since our analysis is done in the framework 
CO ' of Lions [20] and Feireisl et al. [5J , the higher regularity of (p, u) and the uniformly positive lower 

bound of p are not necessary in our analysis and vacuum may be admitted. Indeed, the upper bound 
of the density p plays the essential role in our proof. 
Keywords: Compressible Navier-Stokes equations; decay estimates 
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This paper is devoted to the asymptotic behavior of the solutions to the Navier-Stokes equations of an 
isentropic compressible fluid: 

f p t + div(pu) = 0, 

(1.1) 

\Q . { (pu)t + div(pu ® u) + VP(p) = pAu + (A + p)Vdivu, 

t-h ; 

l/"") | where the density p = p(t,x) and the velocity u = x), u 2 (t, x), ■ ■ ■ ,u N (t,x)) are functions of the 

time t £ (0, oo ) and the spatial coordinate x € Q where fl C M. N , N > 2, is a bounded regular domain. 
| P{p) — ap 1 is the pressure, with a > and 7 > 1 being two positive constants. Since the constant a does 

not play any role in the analysis, we assume henceforth that a — 1. The constants p and A are viscosity 
coefficients satisfying 

p>0, X+^p>0. 

■ We prescribe the initial conditions for the density and momenta: 

p(0) = Po, (pu){0) = m , (1.2) 

together with the no-slip boundary conditions for the velocity: 

u\dn = 0. (1.3) 

Some of the previous works in this direction can be summarized as follows. The first general result on 
weak solutions to the multidimensional isentropic compressible Navier-Stokes equations with large initial 
data was obtained by Lions in |20j . in which he used the renormalization skills introduced by DiPerna 
and Lions in [5] to obtain global weak solutions provided that the specific heat ratio 7 is appropriately 
large, for example, 7 > 3N/(N + 2), TV = 2,3. Later, Feireisl, Novotny and Petzeltovy [5] improved 
Lions's result to the case 7 > If the initial data was assumed to have some symmetric properties, 
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Jiang and Zhang [13l HI] obtained the global weak solutions for any 7 > 1. For the full Navier-Stokes 
equations, Feireisl and Petzeltova established the variational solutions, see [10], for example. 

Concerning the large time behavior of solutions to the initial-boundary value problem (jl.lj) - (jl.3j) . by 
using the weak convergence method, Feireisl and Petzeltova [7 proved the weak solutions to the problem 
(|l.ip - (|1.3[) with a gradient external force VF independent of time t converge to stationary solution (p s , 0) 
in the following sence 



p(t) — > p s strongly in L 7 (0), ess sup / p(t)\u(t)\ 2 dx 

r>t Jq, 



0, as t — >■ 00, 



where the domain Q need not to be bounded and the initial data need not to be close to the equilibrium 
state. If the initial data is close to the equilibrium state, there are many results on the problem of large 
time behavior of global smooth solutions to the compressible Navier-Stokes equations (of heat-conducting 
flow). When there is no external or internal force involved, the LL S global existence and time-decay rate 
of strong solutions are obtained in whole space R 3 first by Matsumura and Nishida [21] [22] and the 
optimal LP (p > 2) decay rate is established by Ponce [25 . The large time decay rate of global solution 
in multi-dimensional half space or exterior domain is also investigated for the compressible Navier-Stokes 
equations by Kagei and Kobayashi [TSUH], Kobayashi and Shibata [T7]> and Kobayashi 18 . Therein, 
the optimal L? time-decay rate in three dimension is established as 

\\(p-p,u)(t)\\ LHm <C(l+t)-i, 

with (p, 0) the constant state, under small initial perturbation in Sobolev space. When additional (exterior 
or internal) potential force is taken into account, the global existence of a strong solution and convergence 
to steady state are investigated by Matsumura and Nishida [24] and many other authors [3] [4] [26] [27j [30] . 
The optimal LP convergence rate in K 3 is established by Duan et al. [5] for the non-isentropic compressible 
flow as 

||(p-p,u,0-0oa)(*)IU»(R 3 ) <C(l + *)~ i(1 ^\ 2<p<6, 

where (p, 0,£?oo) is related to the steady-state solution, under the same smallness assumptions on initial 
perturbation and the external force. If Q is a smooth bounded domain in IR 3 , based on high order energy 
estimates, Matsumura and Nishida 23 proved that, for large times, the solution decays exponentially 
to a unique equilibrium state. However, for the one dimensional case, the smallness assumption on the 
initial data and force can be removed. Indeed, by constructing suitable Lyapunov functionals, the decay 
rate estimates in L 2 norm and H 1 norm are established by Straskraba and Zlotnik [35], and the decay is 
exponential if so the decay rate of the nonstationary part of the mass force is. 

In this paper, under the hypothesis that p is upper bounded, we construct a Lyapunov functional for 
the multidimensional isentropic compressible Navier-Stokes equations (jl.ll) with the aid of the operator 
B introduced by Bogovskii 1 (cf. Lemma 12.51) . Based on this, we show that the weak solutions to 
problem (jl.ip decay exponentially to the equilibrium state in L? norm. The ideas mainly come from [28] , 
however, unlike [28], we do not divide by p on both side of . As a result, the uniformly positive 

lower bound of p is not necessary in our analysis and vacuum may be admitted. Compared with [23], 
our analysis follows the framework of Lions [20] and Feireisl et al. [9], and thus the higher regularity of 
(p, u) is not necessary here. Actually, the upper bound of the density p plays the essential role in our 
proof. Coincidentally, a blow-up criterion for the 3D compressible Navier-Stokes equations was given in 
terms of the upper bound of the density p by Sun et al. [29) . however, their result does not contain the 
case for spatial dimension N > 3. 

Now we give a precise formulation of our result. Let p s be the solution of the following stationary 
problem: 



(1.4) 




Then p s — — — ■ / podx be a positive constant. 
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Formally, the total energy of problem (1.1) can be written as 

E(t)= [ \ P {t)\u(t)\ 2 + -^—p\t)dx, 

satisfying the energy inequality 

— + / p\Vu\ 2 + (\ + p)(divu) 2 dx < 0. (1.5) 
dt Jq 

The definition of weak solutions to the problem (|l.l[) - (ll.3p is given as follows: 

Definition 1.1 (finite energy weak solutions, [20, 9 ). A pair of function (p, u) will be termed a finite 
energy weak solution of the problem (1.1), (1.3) on (0, oo) x Q, if 

• P > 0, p e L% c (0, oo; L-y(Q)), u € Lf oc (0, oo; <' 2 (fi)). 

• The equations (1.1) are satisfied in ^'((0,oo) x Q); moreover, (l.l)i holds in ^'((0,oo) x R N ) 
provided p, u were prolonged to be zero on R^il. 

• The energy inequality (1.4) holds in &'(0, oo). 

• The equality (1.1) i holds in the sense of renormalized solutions, more precisely, the following equa- 
tion 

b(p) t + div(b(p)u) + (b'(p)p - b(p))divu = 
holds in ^'((0, oo) x il) for any b e C X (R) such that 

b'(z) — for all z el large enough, say, \z\ > M, 

where the constant M may vary for different functions b. 

The paper is mainly concerned with the proof of the following theorem. 

Theorem 1.1. Let ft be a bounded Lipschitz domain in l w , N > 2, (p,u) be a finite energy weak 
solution to the problem li._?J) -i TO|) and p s be the solution to the stationary problem |_?.^[ ). In addition, if 
we assume p is upper bounded, i.e., there exists constant p > 0, such that 

p < p, a.e. (t, x) € (0, oo) x Q 

Then 

/ p\u\ 2 + (p- p s ) 2 dx < C(E , p) evp{-C(p, n)t} for a.e. t e (0, oo), (1.6) 



where 



\m Q \ 2 , 1 



E a = / — + -Podx. 

2p 7-1 



Remark 1.1. A natural question is that whether the solution {p, u) stated in Theorem 11.11 exists . Indeed, 
Matsumura and Nishida [23] obtained the existence of global solution to the compressible heat-conductive 
fluid in bounded domain in R 3 provided the initial data is close to the equilibrium state. If vacuum is 
taken into account, Huang, Li and Xin [T^] established the global existence and uniqueness of classical 
solutions to the Cauchy problem for the isentropic compressible Navier-Stokes equations in three spatial 
dimensions with smooth initial data which are of small energy but possibly large oscillations with constant 
state as far field which could be either vacuum or non-vacuum. 

Remark 1.2. We believe that our method can be adapted to the other related models. This is the object 

of our future work. 

Notations: 

1. rj c {-) = -prri(-), where ij is the standard mollifier in IR^. 

2. [/] e = 7 ?e */,forany/eL i 1 oc (K Ar )- 

The rest of the paper is organized as follows: In section 2, we present some preliminary results which 
will be used later. In section 3, we give the proof of Theorem ll.il 
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2 Preliminaries 

Lemma 2.1 ([9]). Let p, u be a solution of (l.l)i in @'((0, oo) x fl) and such that p G L 2 ((0, oo) x 17) 
and it G L 2 (0, oo; [Wo 1 ' 2 ^)]^). 

Then, prolonging p, u to be zero on M. N \Sl, the equation (l.l)i holds in ^'((0,oc) x R N ). 

Lemma 2.2 ([9]). Let (p,u) be a finite energy weak solution of problem (1. !)-(!. 3) on the time interval 
(0,oo). 

Then the total mass m[p(t)} = / p(t)dx is conserved, i.e., 



n 



p{t)dx — I podx, (2-1) 
n Jn 



for all t G (0, oo). 

Lemma 2.3 ([THUD])- Let Cl C R N be a domain and p G L p (fl),u G [W /1 - 9 (ri)] Ar be given functions with 

T7ien /or an?/ compact K a Q, 
ft) 



\\[div(pu)] t - div(\p] e u)\\ L r( K ) < c(K)\\p\\ LP{n) \\u\\ w i, q ( n) (2.2) 

provided e is small enough, where - = - + -. In addition, if fl = M. N , K can be replaced by K. . 
fti) 

[div(pu)] e - d^([p] e w) -> in L e (K) as e -> 0, (2.3) 

where ^ = -^ + ^,ifp<oo and l<0<qifp = oo. In addition, if O = R N and p < oo, if can oe 
replaced by M. N . 

Proof. Since the proof of most of the results in this lemma can be found in |19[ 110) , here we only prove 
(ii) for the case p = oo. To this end, we define G t (p) = [div(pu)] e — div([p] e u). Choosing any open subset 
U such that K C U CC ft, then p G L°°{Q) implies p G IP{U), where P = ^ satisfying i + ± = §. It is 

easy to see that G e (p) — > in L e (f2) as e — ► for any p G Cq°(0). Now choosing a sequence p„ € C§°(U) 
such that p„ — > p in L P (U) as n — >• oo, using the result in (i) with p, p, r and replaced by p — p n , p, 9 
and U, respectively, we have 

\\Ge{p)\\ L C(K) < \\G e (p - Pn)IUe(K) + II Ge(Pn) IU e (if) 

< c(K)\\p - p n \\ LHU) \\u\\ W i, q(u) + \\G e (p n )\\ L e {K) 

< c(K)\\p - p n \\ L p {u) \\u\\ w i, q(n) + \\G e (p n )\\ L e {n) 

— > 0, as e — > 0, n —> oo. 

This completes the proof of Lemma [2T3T □ 

Corollary 2.1. // p G L°°((0, oo) x fi) and u G L 2 (0, oo; [Wq' 2 (CI)] n ) solves (l.l)i m 0'((O, oo) x Q). 
Then for any subset [a, ft] C (0, oo), we have 

dt[p] e + div([p] e u) = r e a.e. on [a, ft] x Q, (2.4) 

where r e — div([p] e u) — [div(pu)] e . Moreover, r e is bounded in L 2 ([a,ft] x Q) uniformly in e and converges 
to strongly in L 2 (a, ft; L 9 (fl)) for all 9 G [1,2). 

Proof. Firstly, by virtue of Lemma 12.11 the equation (l.l)i holds in 2>'({0, oo) x M. ) provided p, u were 
extended to be zero on R w \ri. Then we use the mollifier n e as test functions to deduce (|2.4j) provided e > 
is small enough. It follows from Lemma \2 . 31 immediately that r t is bounded in L 2 ([a, ft] x Q) uniformly 
in e, together with Lebesgue's dominated convergence theorem, we have r e — > in L 2 (a, ft; L e (tt)) for 
9 G [1,2), where we have used the fact that f2 is bounded. □ 
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Lemma 2.4 ([8l [9]). Let be a bounded Lipschitz domain in M. N , N > 2, and p, r £ (l,oo) given 
numbers, f £ / n fdx = 0}. 

Then the problem 

diw — f, v\d£l — 0, (2.5) 
admits a solution operator B : f h- > u enjoying the following properties: 

• B is a linear operator from L p (Tt) into [W ' p (Cl)] N , i.e., 

l|£[/Wi.*(n)||<e(p,Q)||/|U P( n); 

• TTie function v = B[f] solves the problem \2.5\) : 

• If a function f £ L p (£l) can be written in the form f = divg with g £ [L r (fl)] N and g ■ n — on 
dn, where n is the outward pointing unit normal vector field along dtt, then 

\\B[f}\\ Lr{n) <c(p,r,n)\\g\\ Lr{n) . 

Remark 2.1. To our best knowledge, the operator B[-] was first constructed by BogovskiijT]. A complete 
proof of the above mentioned properties may be found in Galdi |11) or Borchers and Sohr 2 . Moreover, 
the operator £>[•] was first used by Feireisl and Petzeltova [8. to show the existence of weak solutions 
{p, u) to the problem (|1.1[) - (|1.3I) with the density p square integrable up to the boundary dfl. 



3 proof of theorem 11.11 

r r ffi _ r T 

Lemma 3.1. Let ro > 0,r > and 7 > 1 be arbitrary fixed constants, f(r) = r I — dh for 

J TO " 

r £ [0, f]. Then there exists positive constants K\ and K2 depending on r$ and f, such that 

Ki(r -r ) 2 < f(r) < K 2 (r - r ) 2 for all r £ [0,f]. (3.1) 

Proof. Let 



l(r) = 



hl - r ldh 



h? 



(r - r ) 2 
It is easy to see that 

l im r / n _ r ° dh 



limj(r) = • -° a " = > 0, 



Using the l'Hospital rule, we obtain 



hl - r l dh+ r lzll 



lim q(r) = lim ■ — ; ; = — rl 2 > 0. 

r^r yy ' r^r 2(f - T ) 2 ° 

Consequently, g(r) is a continuous function on [0,f] with g(r) > 0, and p.lj) follows immediately. □ 

Now we are going to give the proof of Theorem 11.11 First of all, we need to rewrite the energy 
inequality (|1.5j) as a new form. To this end, we choose an arbitrary ip(t) £ 0(0, oo) with ip(t) > 0, then 
the energy inequality (|1.5p is equivalent to 

' f 1 POO r 

ipt / -p\u\ 2 + — — dxdt+ / ip /z|Vw| 2 + (A + /.i)(divu) 2 dxdt < 0. 
o in 2 7-1 Jo in 
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Recalling that p s = — — - j podx is a positive constant and / p(t)dx is independent of t due to Lemma 
l"l Jn Jn 

12.21 we thus have 

1 

7 



It / (-— L 1 pp1 + Pl)dxdt = 0. 
Adding the above two equations, we have 



ipt j (-p|w| 2 + p — ^-dh)dxdt+ I V / (a*|Vm| 2 + (X + p)(divuf)dxdt < 0. (3.2) 



n 



Next, we use the operator B introduced in Lemma [2741 to construct test function of the form 

$(t,x)=V(t)B[[p] e -^], 

where ip is the same as in (|3.2[) and p\ = -r^-r j [p] e dx. Obviously, x) is smooth in x and vanishes near 

c?fi due to the properties of operator £>. Moreover, since p € L°°(0, oo; Q), $ t is in £ 2 (0, oo; [W ' (ty] N ) 
in view of Corollarv l2.ll Consequently, $ could be used as a test function for the equation (|1.1[> 9. Thus, 
we have 

p oo p p oo /> 

V>t / /° u <6[[/>]e — p e s }dxdt + / ?/> / pu$[div([p] e tt)]ctecit 
Jsi Jo Jn 

ip / puB[r t — — - / r e dx]dxdt — / ip pu <g> u : VB[[p] € — p € s ]dxdt 



i> I (P(p) - P(p s ))([p]e - Pl)dxdt + n I p, I Wu: VB[\p] e - p^dxdt 



+(\ + fi) ip divu([p] e - p e s )dxdt = 0, (3.3) 
Jo Jn 

where we have used Corollarv l2.ll Letting e — > 0, we obtain by virtue of Corollary 12.11 and the properties 
of operator B, 

/•OO /> poo P 

V't / puB[p — p s ]dxdt + / ip j puB[div(pu)]dxdt 
Jo Jn Jo Jn 

ip pu ® u : VB[p — p s ]dxdt 
Jo Jn 

/•OO P pOO P 

ip / (P(p) - P( Ps )){p - Ps )dxdt + p / V / Vu : VB[p - p s ]dxdt 
Jo Jn Jo Jn 

+(A + m) / ^ / divit(p - p s )dxdt = 0, (3.4) 



Multiplying (|3.4p by a negative constant — ct with < a <C 1 and summing up the resulting equation 
with the energy inequality (|3.2p . we get 



4>t J (— p|w| 2 + p J — j^-dh- apuB[p- p s ])dxdt 

+ I ip (p\Vu\ 2 + (A + p)(divu) 2 )dxdt - a / ip I puB[div(pu)]dxdt 
Jo Jn Jo Jn 

/•OO P P OO P 

+<t j ip I pu ® u : VS[p — p s ]dxdt + a / ip / {p 1 — pl){p — p s )dxdt 
Jo Jn Jo Jn 

POO P POO P 

—op I ip I : V£>[p — p s ]dxdt — er(A + p) I tp divu(p — p s )dxdt < 0. (3-5) 



Let 



Va = Ji (~2 P ^ 2 + P J — j^-^dh- <rpuB[p- p s ])dx, 



G 



and 

W„ = / (p\Vu\ 2 + (A + fi)(divu) 2 )dx - a j puB[div(pu)]dx 



+a / pu®u: VB[p — p s ]dx + a / (p 1 — p])(p — p s )dx 
Jn Jn 

—up I Vu : VB[p — p s ]dx — a(X + p) I divu(p — p s )dx. 
Jn Jn 



Using the fact p < p and the properties of operator B, we have 



(7 f (7 pc(^l*) f 

-<rpuB[p - p s ])dx\ < — I p\u\ 2 dx-\ — / (p-p s ) 2 dx. (3.6) 

n ^ Jn ^ Jn 



Then it follows from (13.61) and Lemma [37X1 that 



c (a,p,n) / p\u\ 2 + (p- Ps ) 2 dx <V a < Cl {<r,p 7 n) / \u\ 2 + (p- p s ) 2 dx, (3.7) 
Jn Jn 

provided a is small enough. 

On the other hand, by Holder's inequality, we have 

I ~ cr J n puB[div(pu)]dx\ < cr\\pu\\ L 2 (n) \\B[div(pu)} || L 2 (n) 

< crc(n)||p«||| 2(n) < ac(il)p 2 \\u\\ 2 L2{ny 

\a J n pu®u:VB{p- p s ]dx\ < ap( \u\ 2p dx)? (/ | VB[p - p s ] \ q dx) « 

Jn Jn 



<apc{n) [ \\7u\ 2 dx( [ \p-p s \ q dx)' 
Jn Jn 

<<rp 2 c(n) j \Vu\ 2 dx, 
Jn 



{ N if TV > 3 
N ~ 2 ~ ' and - + - = 1. 

2 if JV= 2 p 9 

| -crpj n Vu : Vi3[p-p s ]dx| <j / \Vu\ 2 dx + a 2 p j \VB[p - p s ]\ 2 dx 

4 Jn Jn 

< j / Vw| 2 efe + (T 2 ) uc(^) / \p-p s \ 2 dx, 

| — er(A + /i) / divu (p — p s )dx | < — -^-^ / divit| 2 di£ + cr 2 (A + /i) / \p — p s \ 2 dx, 
Jn 4 Jn Jn 

and it is easy to see that, 

° \ (p 1 ~ Pl)(P ~ Ps)dx >ac \p- p s \ 2 dx. 
Jn Jn 

In view of the above five estimates, we have 

W a >c 2 (a,p,n) [ \u\ 2 + {p- p s fdx, (3.8) 
Jn 

provided a is small enough. 

Therefore, from (|3.7j) and (|3.8|) . we deduce that for appropriate selected cr -C 1, there exists positive 
constant C(p, f2), such that 

C(^,n)^ < (3.9) 

Combining ()3 .9|) with (|3.5[) . we obtain 

oo />oo 

Mt)V«(t)dt + C(p,n) ip(t)V a {t)dt<Q, (3.10) 
o Jo 



for any ip G S>(0, oo) whit tp > 0. 

Let [a, /?] be any compact subset of (0, oo), taking ip(t) = r] e (t — •) in (|3.10|) . we have 

d t [V a ] e + C{p,Q)[V a ] e <0, a.e. i e [a,/3], (3.11) 

provided e is small enough. 
Thus 

WMQ < [VaUs)ex P {-C(p,n)(t - a)} 
for a.e. < s < t < oo, according to (|3.10l) . Recalling that V a (t) G £^ c (0, oo), letting e — > 0, we have 

V a (t) <V a {s)ex V {-C(p,tl)(t-s)} <C(E ,p)exp{-C(p,Sl)t}, (3.12) 

where Eq denotes the initial energy. 

Consequently, (|1.6p follows from (|3.7p and (|3.12[) immediately. This completes the proof of Theo- 
renOJ 

□ 



Acknowledgements 

This work is partially supported by NSFC grant No.10871175, 10931007, 10901137, Zhejiang Provincial 
Natural Science Foundation of China Z6100217, and SRFDP No. 20090101120005. 

References 

[1] M. E. Bogovskii, Solution of some vector analysis problems connected with operators div and grad 
(in Russian), Trudy Sem. S. L. Sobolev 80 (1980), 5-40. 

[2] W. Borchers, H. Sohr, On the equation rotw = g and divu = / with zero boundary conditions, 
Hokkaido Math. J. 19 (1990), 67-87. 

[3] K. Deckelnick, Decay estimates for the compressible Navier-Stokes equations in unbounded domains. 
Math. Z. 209 (1992), 115-130. 

[4] K. Deckelnick, L 2 -decay for the compressible Navier-Stokes equations in unbounded domains, Comm. 
Partial Differ. Equ., 18 (1993), 1445-1476. 

[5] R. J. DiPerna, P.-L. Lions, Ordinary differential equations, transport theory and Sobolev spaces. 
Invent. Math., 98 (1989), 511-547. 

[6] R.-J. Duan, S. Ukai, T. Yang, H.-J. Zhao, Optimal convergence rates for the compressible Navier- 
Stokes equations with potential forces, Math. Models Methods Appl. Sci. 17 (2007), 737-758. 

[7] E. Feireisl, H. Petzeltova, Large-time behavior of solutions to the Navier-Stokes equations of com- 
pressible flow, Arch. Rational Mech. Anal. 150 (1999), 77-96. 

[8] E. Feireisl, H. Petzeltova, On integrability up to the boundary of the weak solutions of the Navier- 
Stokes equations of compressible flow. Comm. Partial Differential Equations, 25 (2000), no. 3-4, 
755-767. 

[9] E. Feireisl, A. Novotny, H. Petzeltova, On the existence of globally defined weak solutions to the 
Navier-Stokes equations, J.Math. Fluid Mech., 3 (2001), 358-392. 

[10] E. Feireisl, Dynamics of Viscous Compressible Fluids, Oxford University Press, Oxford, 2004. 

[11] CP. Galdi, An introduction to the mathematical theory of the Navier-Stokes equations, I, Springer- 
Verlag, New York, 1994. 



8 



[12] X. D., Huang, J. Li, Z. P., Xin, Global Wcll-Poscdness of Classical Solutions with Large Oscil- 
lations and Vacuum to the Three-Dimensional Isentropic Compressible Navier-Stokes Equations, 
arXiv:1004.4749v2 [math-ph] 8 Jul 2010. 

[13] S. Jiang, P. Zhang, On sperically symmetric solutions of the compressible isentropic Navier-Stokes 
equations, Commun. Math. Phys., 215 (2001), 559-581. 

[14] S. Jiang, P. Zhang, Axisymmctric solutions of the 3-D Navier-Stokes equations for compressible 
isentropic flows, J. Math. Pure Appl., 82 (2003), 949-973. 

[15] Y. Kagei, T. Kobayashi, On large time behavior of solutions to the compressible Navier-Stokes 
equations in the half space in M 3 , Arch. Rational Mech. Anal., 165 (2002), 89-159. 

[16] Y. Kagei, T. Kobayashi, Asymptotic behavior of solutions of the compressible Navier-Stokes equa- 
tions on the half space, Arch. Rational Mech. Anal., 177 (2005), 231-330. 

[17] T. Kobayashi, Y.Shibata, Decay estimates of solutions for the equations of motion of compressible 
viscous and heat-conductive gases in an exterior domain in M 3 , Comm. Math. Phys. 200 (1999), 
621-659. 

[18] T. Kobayashi, Some estimates of solutions for the equations of motion of compressible viscous fluid 
in an exterior domain in K 3 , J. Differ. Equ. 184 (2002), 587-619. 

[19] P.-L. Lions, Mathematical Topics in Fluid Mechanics, Vols. 1, Incompressible Models, Oxford Science 
Publication, Oxford, 1996. 

[20] P.-L. Lions, Mathematical Topics in Fluid Mechanics, Vols. 2, Compressible Models, Oxford Science 
Publication, Oxford, 1998. 

[21] A. Matsumura, T. Nishida, The initial value problem for the equation of compressible viscous and 
heat-conductive fluids, Proc. Jpn. Acad. Ser-A., 55 (1979), 337-342. 

[22] A. Matsumura, T. Nishida, The initial value problem for the equation of motion of viscous and 
heat-conductive gases, J. Math. Kyoto. Univ. 20 (1980), 67-104. 

[23] A. Matsumura, T. Nishida, Initial-boundary value problems for the equations of motion of general 
fluids, Computing methods in applied sciences and engineering, V (Versailles, 1981), 389-406, North- 
Holland, Amsterdam, 1982. 

[24] A. Matsumura, T. Nishida, Initial boundary value problems for the equations of motion of compress- 
ible viscous and heat conductive fluids, Comm. Math. Phys., 89 (1983), 445-464. 

[25] G. Ponce, Global existence of small solution to a class of nonlinear evolution equations, Nonlinear 
Anal. 9 (1985), 339-418. 

[26] Y. Shibata, K. Tanaka, On the steady compressible viscous fluid and its stability with respect to 
initial disturbance, J. Math. Soc. Jpn. 55 (2003), 797 826. 

[27] Y. Shibata, K. Tanaka, Rate of convergence of non-stationary flow to the steady flow of compressible 
viscous fluid, Comput. Math. Appl., 53 (2007) 605-623. 

[28] I. Straskraba, A. A. Zlotnik, On a decay rate for ID-viscous compressible barotropic fluid equations, 
J. Evolution Equations, 2 (2002) 69-96. 

[29] Y. Z. Sun, C. Wang, Z. F. Zhang, A Beal-Kato-Majda blow-up criterion for the 3D compressible 
Navier-Stokes equations, J. Math. Pures Appl., in press. 

[30] S. Ukai, T. Yang, H.-J. Zhao, Convergence rate for the compressible Navier-Stokes equations with 
external force, J. Hyperbolic Differ. Equ., 3 (2006), 561-574. 



9 



